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PRESSURE EFFECT OF A SOUND WAVE ON A 
CRYOGENIC FLUID 

SUMMARY 


It has been observed that a sound wave incident upon the surface of a 
liquid in an open container can produce a positive pressure differential between 
the liquid and the ambient atmosphere. This pressure is noticeable only when 
the liquid is close to its boiling point. The proposal has been made to utilize 
such a pressure differential for keeping a cryogenic propellant settled at the 
bottom of its tank during gravity-free flight in space. The report shows that 
this pressure can be explained as an effect of cavitation of the liquid during the 
negative half-cycle of the acoustic pressure wave. 


INTRODUCTION 


When a sound wave encounters a discontinuity in its path of propagation 
such as the boundary between two dissimilar materials, only a part of the 
sound wave will be transmitted into the second medium, while the remaining 
part will be reflected back into the first medium. This reflection of the sound 
energy causes the reflecting surface to experience a pressure which is com- 
monly known as the sound pressure. This sound pressure is a periodic function 
of time according to the periodic character of the sound wave; however, its 
overall, integrated effect on the surface upon which it acts is zero, since the 
effect of the positive half-wave is fully compensated by the effect of the sub- 
sequent negative half-cycle of the pressure wave. It has been demonstrated, 
however, that a positive pressure differential can be observed between a liquid 
in an open container and the atmosphere outside the container, if the liquid is 
near its boiling point and its surface is subjected to an acoustic sound field. 

A typical liquid close to its boiling point may be represented by a cryogen 
under ambient atmospheric pressure. 

This phenomenon of a positive pressure component has been termed 
"acoustical pumping, " and it has been postulated that this positive pressure 
could be utilized to cause residuals of cryogenic propellants to agglomerate 



at a container wall opposite a sound source when the container system is in 
a gravity-free environment. 

This report shows that the phenomenon of "acoustical pumping" can be 
explained by cavitation of the liquid during the negative pressure half-cycle of 
the incident sound wave. 


EFFECT OF SOUND WAVE INCIDENT UPON A FLUID SURFACE 


The apparatus in which the positive pressure component has been observed 
is shown schematically in Figure 1. Before an attempt is made to explain the 



FIGURE 1. SCHEMATIC OF DEMONSTRATION APPARATUS 

observed pressure differential, it is necessary to reiterate the process which 
takes place when a sound wave encounters a discontinuity, such as the interface 
between the air and the liquid in the container. It is known that only a part of 
the sound wave normally incident upon the liquid surface will be transmitted 
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into the liquid. The remaining part of the sound wave will be reflected back 
into the air or vapor on top of the liquid. The intensities of the reflected and 
the transmitted wave are known from the physical properties of the two media; 
the relevant equations are given in Appendix A. These equations have been 
published in the available literature on acoustics [i, 2] , and are given in the 
Appendix to make this report self contained. These equations are valid for the 

idealized conditions that (1) the second medium extends to infinity ox* is 
terminated in such a manner as to resemble acoustically a medium of semi- 
infinite extension; and (2) that the surface of the second medium at the inter- 
face is optically flat and smooth and remains smooth and unaffected by the 
presence of the sound wave. 

In the demonstration apparatus, however, the liquid was only about 5 cm 
(2 in. ) deep as compared 'to the container diameter of about 30 cm ( 12 in. ) , 
or to the wavelength of from 30 to 60 cm ( 12 to 24 in. ) for the employed sound 
frequencies. Also, the surface of the liquid did not remain optically smooth, 
but became rippled and somewhat agitated under the influence of the sound 
vibration. It is evident that the conditions underlying the theory were not fully 
realizable in the demonstration apparatus. The reflection and transmission 
actually taking place in the apparatus must, therefore, to some degree be 
different from the theoretically predicted values. However, for the purpose of 
this report, the theoretical relations as given by the equations in Appendix A 
are accepted also for conditions as they prevailed in the demonstration appara- 
tus. As a consequence, the incident wave should be almost totally reflected 
from the liquid surface, and the pressure experienced by the surface of the 
liquid is expected to be twice the pressure in the incident wave. 

If a liquid is subjected to a pressure, it will transmit this pressure as an 
omnidirectional stress through its entire volume. This is basically the inter- 
pretation of equation (A-9) when viewed from the field of fluid mechanics. 

It is also known that a liquid has an unlimited capacity for supporting compres- 
sive stresses (positive gage or absolute pressure) , but that it is limited in its 
capability to transmit a tensile stress (or negative gage pressure) . This limit 
is established by the vapor pressure of the liquid. If the tensile stress in the 
liquid tends to be lower than the vapor pressure, the liquid will tear or 
separate, and form cavities. These cavities are filled with saturated vapor of 
the liquid. This formation of vapor-filled cavities is known as cavitation. 

The vapor pressure is a function of temperature. The relation 
p^ = f(T) is shown for water in Figure 2. It shows that at 372. 25® K (99. i®C) 

water has a vapor pressure of 1 kg/cm 2 or one atmosphere absolute. At 
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Pressure 
N/m 2 x 10 -5 

I 

1. 961 3300 

1.765 1970 
1.569 0640 
1.372 9310 
1.176 7980 


0. 980 6650 
0. 882 5985 
0. 784 5320 
0.686 4655 
Vapor pressure of water 

k = m > I 


413 433 


FIGURE 2. VAPOR PRESSURE OF WATER AS A 
FUNCTION OF TEMPERATURE 


453 
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354. 15* K ( 81*C) the vapor pressure is only 0. 5 kg/cm 2 or one-half atmosphere. 
Therefore, water at 372. 25* K (99. 1*C) cannot support any tensile stress (or 
negative gage pressure below one atmosphere ambient) but it can sustain a 
tensile stress or negative gage pressure of -0. 5 kg/cm 2 if its temperature is 
lowered to 354. 15* K (8i*C). If the attempt were made to lower this pressure 
to a more negative value below -0. 5 kg/cm 2 , the liquid would again cavitate, 
that is, form cavities filled with saturated vapor at 354. 15*K (81°C); the 
pressure in the liquid would not change in spite of the added vapor volume. 

If a periodically varying pressure is superimposed to the ambient 
pressure which prevails on the surface of a liquid, a total absolute pressure 
on the liquid's surface results as a function of time as shown in Figure 3(a) . 

If the temperature of the liquid is low, the vapor pressure pertaining to this 
low temperature will also be well below the absolute pressure acting on the 
surface. 

If the temperature of the liquid is now raised, the vapor pressure of the 
liquid will also rise until it falls within the negative half-wave of the super- 
imposed pressure fluctuation. This is illustrated in Figure 3(b). Since the 
liquid cannot sustain a negative pressure below its vapor pressure, cavitation 
must now occur during that part of the negative pressure cycle which falls 
below the value of the vapor pressure. This is indicated by the shaded part in 
Figure 3(b) . It is also to be noted that the negative pressure amplitude below 
the vapor pressure within the shaded part of the negative half-wave is not 
transmitted into the liquid; therefore, a positive pressure component above 
ambient results from the integral Jp • dt taken over any number of complete 
cycles. 

In the extreme case illustrated in Figure 3( c) , the temperature of the 
liquid has been raised to the boiling point. In this case, only the positive 
pressure can be transmitted through the liquid; the negative pressure below 
ambient causes only cavitation or a release of additional vapor in excess of the 
boil-off rate caused by external heating. In this case, the positive differential 
pressure between the liquid in the container and the ambient pressure outside 
the container will attain its maximum value. 

Cavitation or vapor formation in the examples cited above takes place 
only at the surface of the liquid. This is caused by the gravitational field, which 
causes the pressure in the liquid to increase toward the bottom of the container. 
Any vapor release caused by the reduction of the pressure on the liquid's sur- 
face must, therefore, take place at the surface where the pressure in the liquid 
is closest to the vapor pressure. 
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FIGURE 3. SOUND PRESSURE WAVE SUPERIMPOSED TO ATMOSPHERIC 
PRESSURE IN RELATION TO VAPOR PRESSURE OF COLD LIQUID, 
HOT LIQUID, AND BOILING LIQUID 



ANALYSIS OF EFFECT OF CAVITATION ON PERIODIC 
PRESSURE VARIATION 


The periodic pressure variation of a perfect sound wave can be expressed 
mathematically by a trigonometric function as a pure sine wave. This is not 
possible for the pressure oscillation acting on the surface of the liquid if the 
negative half-wave is only partly transmitted and effective because of the 
cavitating liquid. However, the clipped-off pressure wave shown in Figure 3(a) 
and (b) can be represented by a Fourier series expansion, which also allows 
analysis of the particular characteristics and side effects of such a disturbed 
pressure fluctuation. This has been done in Appendix B for the pressure cycle 
shown in Figure 4. 



FIGURE 4. PERIODIC PRESSURE FUNCTION LIMITED BY CAVITATION 

The analysis shows that this particular wave shape produces a constant, 
positive pressure component, the magnitude of which depends only on the 
reduction of the negative pressure peak by the vapor pressure. It also shows 
that an infinitely large number of higher harmonics is now introduced. These 
higher harmonics may produce special effects in the vibration pattern of the 
liquid and the container, and may actually affect the positive pressure component 
beyond that magnitude which is given by the constant term. Additional cavitation 
may be caused by the higher frequency harmonics because of their high accelera- 
tion levels. These effects, however, are disregarded here. 
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NUMER I CAL ANALYS I S OF ACOUSTI CAL PRESSURE 
CAUSED BY CAVITATION 


A loudspeaker may be used as a source of acoustical energy in the test 
arrangement illustrated in Figure i. The electrical energy, W, supplied to 
the terminals of the speaker is not converted fully into acoustical energy 
because of some losses inherent in the energy conversion. These losses are 

1. Copper losses (electrical energy converted into 

heat in the speaker coil) 

2. Hysteresis and induction losses (caused by the rapidly changing 

magnetic field) 

3. Mechanical losses (caused by the deformation of the 

speaker membrane and its internal 
material damping) 

4. Losses from the non-isentropic boundary layer attached to the 
membrane surface. 

There might be additional losses. The actual efficiency of the sound source 
used in a test arrangement will have to be determined by calibration. 

The sound intensity generated by the loudspeaker is then 


I 

s 



W 

D 2 7T 

4 


W/m 2 


where D = diameter of speaker membrane (m) 

W= electrical energy fed into the speaker (W) 

tj = conversion efficiency ( -) 

c 

The sound energy of a plane, free progressive wave is 


I 


W 


Po 


Po ' c 


W/m 2 


U) 


( 2 ) 
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where p 0 = rms value of the sound pressure 
p 0 • c = characteristic impedance 


(N/m 2 ) 
(N-sec/m 2 ) . 


The sound energy generated by the speaker will be radiated primarily 
into the interior of the container, its magnitude depending on the directivity 
characteristics of the speaker. Part of the sound energy remains on the out- 
side of the system as indicated in Figure i. The ratio of these two energies 
may be denoted by e. The sound energy inside the container is then given by 


T _ W 

in € % D 2 7T 
4 

and the sound intensity outside the container is 


(3) 


I = ( 1 - e) 
out 


W 

D 2 7T 


(4) 


Without sound attenuation, the rms value of the sound pressure inside the 
cylindrical container on the surface of the liquid is obtained from equations 
(2) and (3) as 


Po = 


4e- n 
c 

7T 


-,1/2 


w 

(PoC) 


(5) 


Introducing the numerical value of the impedance for air at 295. 15° K 
(22° C) and 751 mm Hg pressure as Po • c = 407 N-sec/m 2 results in 


€ • 



22 . 7 - 


( 6 ) 


The diameter, D, initially taken as the diameter of the speaker membrane, 
is now to be taken as the inside diameter of the container; however, the two 
diameters should not differ widely from one another. 


It is shown in Appendix A that because of the almost complete reflection 
of the incident sound wave from the surface of the liquid, the instantaneous 
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sound pressure p(t) is twice the pressure of the incident wave. However, the 
pressure amplitude of the incident wave is n/ 2" times the rms pressure p 0 . 
Thus, the amplitude of the pressure acting on the surface of the liquid is 


A = 2 • nT 2 • p 0 = 2. 825 p 0 
Combining equations (6) and (7) finally results in 


1/2 64 D 


(7) 


( 8 ) 


e • n • W 
'c 

Equation ( 8) states that 1 W of electrical energy supplied to a loudspeaker 

of efficiency r? and a diameter of 1 m would produce a pressure amplitude of 
c 

64 N/m 2 . Since a barometric pressure of 751 mm Hg equals 10 5 N/m 2 , the 
pressure amplitude expected at the surface of the liquid, expressed in mm Hg, 
is 

1/2 

mm Hg. (9) 


A Hg“ U 'V • °- 482 \- 


Expressed in head of water, the pressure amplitude is 

1/2 


W 


A W = U * V ’ 6 * 55 D 


mm water. 


( 10 ) 


This value of A, if introduced into equation (B-10) , yields the positive 
pressure component which can be expected as the "acoustical pumping" pressure 
for a given electrical input into the speaker 


p+ = 


0. 482 | * € * W 


1/2 


D 


£sin (/>-(/) cos mm Hg. (11) 


The value of cos <p is found from Figure 4 as 


P - p 
a v 

cos </> = ; 0 < (p - p ) < A 

A *a *v 


( 12 ) 


Using equation ( 8) , the final value is 

D 


COS 4 > = 


64 • € • i) *W 

'c 


T/2 (p a- p v> 


(13) 
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CONCLUSIONS AND RECOMMENDATIONS 


The explanation that the observed positive pressure in the liquid may be 
caused by cavitation seems to be supported by the analysis of the pressure wave 
with clipped-off negative pressure amplitude. It would be of great interest to 
verify the numerical values of pressure obtained from 1 W of electrical energy 
by performing some well controlled but simple tests. The frequency of the 
imposed sound wave should be in the audible range. The pressure observed in 
such a test could then be considered an indication of the validity of the theory 
developed in this report. 

It would also be of great interest to determine the effect of the sound 
frequency of much shorter wavelengths. According to the results shown in 
Appendix B, the frequency should have no influence on the magnitude of the 
pressure. However, it is known that ultrahigh frequency sound has some 
directional characteristics which may cause a deviation from the results shown 
in Appendix B. The extension of the sound frequency into the ultrasonic range 
would, therefore, be very desirable. 

Generally, the omnidirectional pressure in the liquid, as generated by 
the sound wave, will not be sufficient in itself to keep the liquid settled at a 
predetermined location in the tank. The effect of the acoustical pressure on the 
liquid is the same as a commonly performed tank pressurization from an exter- 
nal pressure source. It can not cause the liquid to stay agglomorated at a 
definite place in the tank under a gravity-free environment. To accomplish this 
requires a pressure gradient within the fluid, similar to the static pressure 
increase with depth in the liquid when subjected to gravity, or similar to the 
pressure distribution within a rotating fluid due to centrifugal forces. Such a 
pressure gradient cannot be created by the acoustic pressure investigated in 
this report. However, sound waves with ultrahigh frequencies may produce a 
pseudo pressure gradient due to their directional characteristics. For this 
reason, tests in the ultrasonic frequency range appear to be very desirable. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, April 25, 1967 
103-19-05-00-62 
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APPENDIX A 


REFLECTION OF A PLANE WAVE FROM A 
PLANE BOUNDARY BETWEEN TWO MEDIA 


The plane boundary between media I and II establishes the discontinuity 
at x = 0 in Figure A-i. If p indicates pressure, u the particle velocity, and R 
the characteristic impedance, a simplified analysis [1,2] can be performed 
using the wave equation for sinusoidal waves: 


I ii 



-x x=o +x 

FIGURE A-l. REFLECTION OF PLANE WAVE AT A PLANE BOUNDARY 


p^ = Aj • sin (cot - kjx) = u. • Ri 

(A-l) 

p^ = Bj sin (cot - kjx) = - u^ ' Ri 

( A-2) 

P t = A 2 sin (cot - k 2 x) = u t • R 2 , 

(A-3) 


where 

~ = p • c = R , (A-4) 

Aj, Bj, and A 2 are the pressure amplitudes, 
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k t and k 2 are the wave numbers defined as — , 

A 


and x is the positive direction of the incident wave. 

The subscripts i, r, and t refer to the incident, reflected, and transmitted 
soundwave, respectively. 


At the boundary, the pressure and the velocity must be continuous, thus 


p t = p i + p r 

(A-5) 

U t = U i + U r 

(A-6) 

If equations (A-4) and (A-5) are introduced into equations (A-i) , 
(A- 3) , the following relations are obtained: 

(A-2) , and 

A 2 = A* + Bj 

(A-7) 

Ri * A 2 = R 2 (Ai - Bi) 

(A-8) 

A 2 2 R 2 

p. Aj Rj + R 2 

(A-9) 

P r ^2 ~ 

P j i R i + R 2 

( A-10) 

For the particle velocities follows 


u A 2 Rj 2 R^ 

u. A| R 2 R a + R 2 

(A-li) 

u r B i R i " r 2 

u. Aj Rj + R 2 

(A-12) 


The ratios given in equations (A- 9) through (A-12) are the amplitude ratios of 
particle pressure and particle velocity, respectively. 
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The transition of sound energy from one medium into the other is defined 
as the ratio of the intensity of the transmitted wave to that of the incident wave. 
Intensity is defined as 


from which follows 

I t 4R 1 R 2 

a t = r = rTpT = <Ri + r 2 ) 2 

and the ratio of the reflected intensities is obtained as 

i r P 2 r • Ri <R 2 - Ri) 2 

°r = i~ "rTp 2 = (Ri + R 2 ) 2 


(A-13) 


(A-14) 


(A-15) 


Equations (A-9) through (A-15) show that if R 2 and R 2 are equal, there is no 
sound reflection; the entire sound energy is transmitted into the second medium. 

If both materials are very dissimilar, such that R 2 §> Rj (as is the 
case if medium I is air or saturated steam, and medium II is water) , then the 
numerator in equation (A-12) will be negative. This means that the particle 
velocity changes in phase on reflection; that is, the reflected wave is opposite 
in phase to the incident wave. The phase of the pressure, however, remains 
unchanged. If R 2 approaches infinity, as in the case of an absolutely rigid sur- 
face , the phase of the particle velocity of the reflected wave is opposite to that 
of the incident wave by ir , and its amplitude is equal to the amplitude of the 
incident wave; the resultant particle velocity at the boundary is therefore zero. 
However, since the phase of the pressure does not change, it follows that the 
pressure at the boundary is doubled. Thus, the reflection results in a standing 
wave which has a node at the wall, and a pressure amplitude at the wall of twice 
the amplitude of the incident wave. 
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APPENDIX B 


ANALYSIS OF PRESSURE OSCILLATION 
LIMITED BY CAVITATION 


A periodic function such as shown in Figure 4 can be expressed by a 
Fourier series expansion of the form 


a° n = °° 

f(x) = — + ( a n COS + \ sin ( Qx )) ’ (B-i) 

n = 1 


where 


1 


a 0 7r 


+7T 

J f(x) dx 

-7T 


a 

n 


1 

IT 


+7T 

/ f(x) cos nx dx , 

“7T 


and 


b 

n 


1 

7T 


+7T 

/ fix) sin nx dx . 
-7 r 


( B-2) 


( B-3) 


( B-4) 


The undisturbed pressure wave is of the form -A* cos cut. Cavitation may take 
place between the values of cut! = -<p and cut 2 = +<P at a vapor pressure p y which 

limits the negative pressure half-wave to a magnitude of -A* cos cp . The 
Fourier coefficients are found as follows: 

1 +7T 

a o = 7 r / f (cot) d(cut) 

-IT 


1 


~ 7 r 



- A cos cut d(cut) 




- A cos (p d(cut) 


+ 


1 

n 


7T 


I 

* 


A cos cut d(cut) 
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or finally 



[ 


sin cot 


-4> 

+ COS 4 > cot 

—IT 


+<f> 

+ sin cot 

-<P 



a o= ~ 



( B-5) 


Likewise , 


a 

n 


7T j 

/ f(cot) • cos n cot • d(cot) =|}f J -Acos cot* cos n cot* d(cot) 


-7T 


1 +0 7T 

+ 7T / - A cos (p cos n cot d(cot) + J - A • coscot - cos n cot - d(cot) , 

~<f> + 4> 

which yields 

a = F f _ sStlaiilAl . <b-6> 

n " |_n * (n - i) n • (n + 1) J 


The integrals leading to equation (B-6) are not applicable to the case of 
n = 1; the value of aj (for n = 1) has to be determined from the original Fourier 
equation for a by introducing n = 1 before the integration is performed. This 
yields 


j ~ 4 > 1 +0 

a i = ~ f - A cos 2 cot • d(cot) + — f - A cos <p • cos cot • d(cot) 

~<t> 


— 7T 


or 


+ ^ J - A cos 2 cot- d(cot) , 

+<P 

a^ = - ~ £ 7r - <£ + cos <p • sin (p J 


(B-7) 
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Combining equations (B-5) , (B-6) , and (B-7) , the Fourier series 
expansion for the pressure function shown in Figure 4 reads 


P(t) = 


_ A 

TT 


(9 cos 4> ~ sin i p) + (tt - 9 + cos <j> sin 9) • cos cot + 


Z i sin (n - 1) 9 sin (n + 1) 9 
„ \n (n - 1) 


n (n + 1) 


cos 


( B-8) 


n = oo 

2 

n = 2 

Equation (B-8) shows that for 9 = 0 (no cavitation) , the pressure is given by 


n • cot J . 


p(t) , = - A • cos cot , 

9 = 0 


(B-9) 


which is the original, undisturbed sound pressure oscillation. 

In case of cavitation, 9^0, there is a positive, constant pressure 
component independent of the sound frequency, but dependent only on the 
maximum pressure amplitude A and the cavitation angle <p. This constant 
pressure is given by 


p+ = y (cf> cos 9 - sin 9) 


( B-10) 


The magnitude of this pressure component, p+/A is given in Table B-I and is 
illustrated in Figure B-l. 


TABLE B-I. MAGNITUDE OF POSITIVE PRESSURE 
p+/A = ^ (sin 9 - $ cos 9) 


0 

P+/A 

0 

0. 00 000 

tt/12 

0. 00 623 

2 7T/12 

0. 04 707 

3 ir/12 

0. 15 203 

4 tt/12 

0. 34 253 

5 7r/ 12 

0. 62 701 

6 tt/12 

1. 00 000 
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FIGURE B-l. MAGNITUDE OF POSITIVE PRESSURE COMPONENT 

p+/A = ^ ( sin <f> - <p cos <p) 


In case of cavitation, higher harmonics of the pressure wave are also 
created with pressure amplitudes as given by the term 


P 


n 



sin (n - 1) 0 
n (n - 1) 


sin (n + 1) 0 
n • (n + i) 


cos n • cut . 


The magnitude of the pressure amplitudes of the first five higher harmonics 
(n = 2, 3, 4, 5) are given in Table B-II for the values of (p = n/6, ir/3, 7r/2, 
where the case of 4> = tt/2 corresponds to the case of the boiling liquid. 

TABLE B-2, PRESSURE AMPLITUDES OF HIGHER HARMONICS 


"V n 

9 

2 

3 

4 

5 

7t/6 

- . 02 6526 

- . 02 2972 

- . 01 8568 

- . 01 3783 

7r/3 

- . 13 7833 

- . 06 8917 

- . 01 3783 

+ . 01 3783 

7t/2 

- .21 2207 

0. 0 

- . 04 2441 

0. 0 
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Figure B-2 shows how closely the true wave shape of the three cases for 
4> = Jf/6, 7t/ 3, and n/2 is represented by equation (B-8) , when only four higher 
harmonics (n = 5) are used. 



FIGURE B-2. COMPARISON OF EXACT WAVE SHAPE WITH 
APPROXIMATION BY FOURIER SERIES EXPANSION 

(for n = 5) 
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